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The key process in the climate change fingerprinting is to attribute the 
large spatiotemporal averaged spectral variation to different climate variables. 
This attribution can be expressed as a multivariable linear regression problem: 

AR=KAX+e (la) 


AR^^KtjAXj + ej (lb) 



dx j ’ 


i j l,2,...,n x 


(lc) 


AR Reflectance change signal (spectral difference between two climate states) 
K Kernel matrix (fingerprints) 

AX Climate variable changes to be retrieved 
e Errors or residuals that cannot be explained by fingerprints 

K cannot be measured directly, has to be generated from RT modeling. 

R either from observation or model simulation. 


Without the error consideration, the solution of Equation (1) is simply 


XX = (K T K)~' K T AR 


This is based on the ordinary least squares estimation (LSE). A small error in K 
or in AR could produce large error in the solution. 

A common approach used to account the error is the so-called optimal 
detection: 


AX - (K T E~' K)~' K T E~' AR 


where E is the covariance matrix of e in Equation (1). 




The formulation for fingerprinting retrieval here is similar to those used in 
many conventional retrieval methods applied to instantaneous satellite data. 
However, the fingerprinting retrieval differs from the traditional remote sensing 
retrieval in that 

1) It uses the average-then-retrieval approach and thus is associated with the 
averaged quantities over large spatiotemporal scales instead of the local or 
instantaneous values. 

2) It is for the difference (AR and AX) between two mean climate states instead 
of the absolute value (R and X). 


Here we use model-simulated data to test the concept of fingerprinting 
attribution. Model input parameters: 

S 10 years (2001-2010) of CERES SSF data from NASA Terra satellite ; 

S aerosol and cloud properties (optical depth, particle size, phase, and 
height) retrieved from MODIS; 

S column water vapor and surface wind data from GEOS5-MERRA 
reanalysis; 

S ozone data from SMOBA; 

S ocean chlorophyll concentration from SeaWiFS. 

Using these input parameters in combination with COART-MODTRAN, we 
generate the spectral kernels and a time series (10 years) of monthly mean 
reflectance spectra (320-2300nm, 4nm resolution) to retrieve the interannual 
changes in the 11 relevant climate parameters: 

PI/I/, AOD, 03 

and 

r, Fc, Ht, Re for water and ice clouds, respectively . 
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An example of the monthly and global mean reflectance anomalies (AR) in 
the four months. In each panel, each color is for a different year. 


AR is indeed small: typically less than ±3% of the mean reflectance 



PW 


ice 


- wat 



CUO 

CD 

CD 

“D 


"■P 

CD 


0.5 


1.0 1.5 

Ice Cloud r 


2.0 


0.5 


1.0 1.5 

Water Cloud r 


2.0 




10.000 , 
- 0.002 
- 0.003 
- 0.005 
- 0.006 
- 0.008 
1 - 0.010 
- 0.01 1 
- 0.013 
I - 0.014 
- 0.016 
- 0.018 
- 0.019 


8:8ee 

0.056 
0.054 
0.052 
0.050 
0.048 
0.046 
0.044 
0 . 042 ' 
0.040 
0 . 038 ' 
0.036 
0.034 
0.032 
0.030 
0.028 
0.026 
0.024 
0.022 
0.020 
0.018 
0.01 6 
0.014 
0.012 
0.010 
0.008 
0.006 
0.004 

8:888 


0.5 1.0 1.5 2.0 

Wavelength (pim) 


An example of solar spectral reflectance kernels 


Reflectance Kernel (1/c 
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Comparison between the 
fingerprinting retrieval and the 
observational truth for the 
monthly global mean anomaly. 

Each panel is for a different 
climate parameter (11 total). 

The nonlinear error in kernels is 
not considered in the retrieval. 
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Same retrieval as above, but 
the nonlinear error in kernels is 
considered. 


Two nonlinear error sources: 


1). nonlinear radiative response : 


f; 


x 0 +Ax 


K{x)dx * K(x 0 ) Ax 


2). radiative interactions: 
ARtlj ARj = I j KAXj 
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Without nonlinear errors considered 
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With nonlinear errors considered 



Latitude 


With the averaged kernels and reflectance in the 5 large latitude regions, 
the fingerprinting retrieval as above is tested in these large regions and the 
results are similar as above. 


Five latitude regions 
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RMS comparison in all regions 
fro all 11 parameters. 

Each panel represents a 
different region. 

The RMS is normalized by the 
variation a for each parameter. 
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Without nonlinear error considered (LSE) 
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With nonlinear error considered (optimal detection) 
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The kernel spectral shapes 
for the cloud fraction and 

the cloud G are similar in 
many cases! 

The upper panel is for water 
cloud and the lower panel is 
for ice cloud. 

This example is for a zonal 
mean in the tropic (0N-10N). 
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A comparison of the cloud i 
retrievals between two 
cases: 

1) the cloud fraction is 
unknown and is also 
retrieved (the left 
panels); 

2) the cloud fraction is 
known a priori in the 
retrieval (the right 
panels). 

The upper panel is for water 
cloud and the lower panel is 
for ice cloud. This example is 
for the mean data between 
60N to 60S in latitude. 
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Comparison of the RMS 
errors in the retrieved 
cloud x anomalies in the 
five climate regions 
between known and 
unknown cloud fraction. 
The upper panel is for 
water cloud and the lower 
panel is for ice cloud. 


The improvement is general in all regions! 

Suggesting a separate measure of cloud amount from CLARREO to 
unscramble cloud fraction and optical depth. 


AOD Retrieval 
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A retrieval comparison 
for the three 
atmospheric parameters 
(PW, AOD and ozone) 
between all sky and 
clear sky retrievals. Each 
row is for a different 
parameter. 
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Comparison of the retrieval 
RMS errors in the three 
parameters between clear sky 
and all sky retrievals in the 
five climate regions. 

This example is for the mean 
data between 60N to 60S in 
latitude. 


Conclusion 


S Using ten years of satellite data, we tested the concept of using fingerprinting 
approach for climate change detection/attribution. 

S Comparing the fingerprinting retrieval to the observational truth, the RMS 
differences are less than 2 a of the variance for all variables in all regions. A large 
error usually corresponds to the variables with large nonlinear radiative response, 
such as cloud x and ice particle size. 

S Using the optimal detection method to take into account the nonlinear radiative 
error in the kernels, the retrieval accuracy is significantly higher, so that the RMS 
errors are reduced to less than 1 o of the variance, indicating the profound impact 
of the nonlinear error on fingerprinting retrieval. 

S Another important finding is that if the cloud fraction is known a priori, the 
retrieval accuracy in cloud x would be improved substantially and generally; in 
addition, a better retrieval for the water vapor amount and aerosol optical depth 
can be achieved from the clear sky data only. 

S This finding suggests a separate measure of cloud fraction by CLARREO to 
unscramble cloud fraction and optical depth variations. 


S Ozone retrieval has the highest accuracy among all variables in all cases. 

S The test results demonstrate that the concept of climate change fingerprinting 
based on the reflected solar benchmark spectra is viable. 



